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Abstract. Meynet & Arnould (1993a) have suggested 
that Wolf-Rayet (WR) stars could significantly conta- 
minate the Galaxy with ^^F. In their scenario, ^^F is syn- 
thesized at the beginning of the He-burning phase from 
the ^''N left over by the previous CNO-burning core, and 
is ejected in the interstellar medium when the star en- 
ters its WC phase. Recourse to CNO seeds makes the ^^F 
yields metallicity-dependent. 

These yields are calculated on grounds of detailed stel- 
lar evolutionary sequences for an extended range of initial 
masses (from 25 to 120 Mq) and metallicities {Z = 0.008, 
0.02 and 0.04). The adopted mass loss rate prescription 
enables to account for the observed variations of WR pop- 
ulations in different environments. 

The ^^F abundance in the WR winds of 60 M© model 
stars is found to be about 10 to 70 times higher than its 
initial value, depending on the metallicity. This predic- 
tion is used in conjunction with a very simple model for 
the chemical evolution of the Galaxy to predict that WR 
stars could be significant (dominant ?) contributors to the 
solar system fluorine content. We also briefly discuss the 
implications of our model on the possible detection of flu- 
orine at high redshift. 

Key words: Nuclear reactions, nucleosynthesis, abun- 
dances - Stars: Wolf-Rayet 



1. Introduction 

For long, the solar system has been the only location of 
the Galaxy with a known fluorine (^^F) abundance. At the 
same time, the production site(s) of this element has been 
a major nucleosynthetic puzzle, even if F is the least abun- 
dant (mass fraction of 4 10~^, following Grevesse & Sauval 
1998) of the elements ranging from carbon to calcium. 

These last years, the situation has changed quite dra- 
matically, both observationally and theoretically. Fluorine 
overabundances (with respect to solar) in MS, S and C 



Send offprint requests to: G. Meynet 



stars have been reported (Jorissen et al. 1992), and cor- 
relate in particular with s-process enrichments. These ob- 
servations demonstrate that thermally pulsating Asymp- 
totic Giant Branch (AGB) stars are fluorine producers, 
as predicted by Goriely et al. (1989), and confirmed by 
calculations conducted in the framework of detailed AGB 
models (Forestini et al 1992, Mowlavi et al. 1996, 1998). It 
remains of course to determine the exact level of the con- 
tribution of these (mass losing) stars to the solar system 
and galactic F content. 

Iir direct relation with this question, various calcula- 
tions have been made in order to estimate the ^^F yields 
from massive stars. The neutrino process operating dur- 
ing supernova explosions has been envisioned as a possible 
producer of primary ^^F (e.g. Woosley & Weaver 1995). 
On the other hand, Meynet & Arnould (1993a) have in- 
vestigated on grounds of detailed stellar models the sug- 
gestion (Goriely et al. 1989) that the hydrostatically burn- 
ing He-shell can synthesize ^^F of secondary nature. They 
find that the level of production is relatively modest in 
M < 20 M0. In contrast, they show that stars which are 
massive enough to become Wolf-Rayet (WR) stars can 
eject through their winds substantial amounts of fiuorine 
synthesized in the core at the beginning of the He-burning 
phase. 

In the present work, we revisit the question of the 
galactic contribution of WR stars to ^^F with the help of 
new stellar models that better account for many important 
observable properties of WR stars. In addition, we extend 
the range of masses and metallicities considered in our 
previous study. The broadening of the explored metallic- 
ity range may take some additional importance in relation 
with the recent claim by Timmes et al. (1997) that "posi- 
tive detection of any fluorine at a sufficiently large redshift 
(z > 1.5^ would suggest strongly a positive detection of the 
neutrino process operating in massive stars'^ . The possibil- 
ity of a significant thermonuclear production of ^^F by WR 
stars of different metallicities might blur this picture, and 
might at least imply the necessity of establishing observa- 
tionally the primary or secondary nature of the detected 
fiuorine, if any. 



2 



G. Meynet & M. Arnould: Synthesis of ^''F in Wolf-Rayet stars 



The physical ingredients of the models arc discussed 
in Sect. 2. Section 3 presents our predicted yields from 
individual WR stars, while Sect. 4 gives a rough estimate 
of the contribution of WR stars to the galactic ""^^F content. 
Some conclusions are drawn in Sect. 5. 

2. The physical ingredients of the stellar models 

The evolutionary models are computed with the same 
physical ingredients as in Meynet et al. (1994). However, 
the adopted nuclear reaction network is extended, espe- 
cially in order to include the reactions involved in the 
production and destruction of ^^F (see below). 

The differences with respect to the computations of 
Meynet & Arnould (1993a) arc twofold: 

1) A more extended range of initial masses (from 25 to 
120 Mq) and metallicities {Z = 0.008, 0.02 and 0.04) is 
explored; 

2) The present grid of models is computed with the mass 
loss rates adopted by Meynet et al. (1994), which are twice 
as large as the values of M recommended by de Jager et 
al. (1988) and Conti (1988) for the pre-WR and WNL 
phases. This mass loss rate prescription enables to account 
for the observed variations of WR populations in different 
environments (Maeder & Meynet 1994). 

The metallicity dependence of the mass loss rates dur- 
ing the prc-WR phases is adopted from previous works 
(e.g. Meynet et al. 1994). More specifically, M scales with 
metallicity Z according to Mz/Mq = (Z/Z©)"-^, where 
Zq is the solar metallicity. This scaling is deduced from 
stellar wind models (cf. Kudritzki et al. 1987, 1991). 

Let us finally add that the models are computed with a 
moderate core overshooting (d/Hp = 0.20, where d is the 
overshooting distance and Hp the pressure scale height at 
the boundary of the classical core). 

2.1. The thermonuclear ^^F production and destruction 

paths 

The CNO mode of H-burning is responsible for the pro- 
duction and destruction of ^^F through the reaction chain 

"N(p, 7)150(/3+)15N(p, 7)160(p, 7)"F, 
i^F(/3+)i70(p,7)^'F(/?+)i8o(p,7)i9F(p,a)i6o. 

The adopted ^^F{p, a)^^0 rate is the geometrical mean 
of the lower and upper limits to that rate proposed by 
Kious (1990). 

Fluorine can also be produced and destroyed dur- 
ing He-burning through the chains (see also Meynet & 
Arnould 1993a) 



{(3+y'0{p,ar^N{a,^y'F 



/ \ 
"iV(a,7)i«^^ ^ (n,p)i«0(p,a)i5iV(a,7)i3f^(a,p)22iVe. 

\ / 

The synthesis of ^^F thus requires the availability of neu- 
trons and protons. They arc mainly produced by the re- 
actions i-'^C (a , n) and i4N(n,p)i4c. 

The first chain of transformation of ^''N into ^^F men- 
tioned above is by far the most important in the condi- 
tions of relevance in this work, where the /3+-decay lifetime 
T^i^'^Y) of i*F is much shorter than its lifetime r„,a(i*F) 
or T„,p(i^F) against (n, a) or (n,p) reactions. For example, 
T/3 is a few hours only at the center of a 60 M© model at 
the beginning of core He-burning, while the corresponding 
T„,a(^^F) and r„,p(i^F) amount to about 1 400 and 18 000 
years, respectively. 

The NACRE compilation of reaction rates (Angulo et 
al. 1999) was not available yet at the time of completion 
of the calculations reported here. This is why most of the 
necessary nuclear data are taken from Caughlan & Fowler 
(1988). There are some exceptions to this rule. In par- 
ticular, the I'^C a-capture rate is taken from Descouve- 
mont (1987), whose theoretical prediction of an increase 
of the astrophysical S-factor at low energies is confirmed 
experimentally (see NACRE). The I'^N (n , p) ^''C rate is 
taken from Brehm et al. (1988). It is a factor of two lower 
than the one proposed by Koehler and O'Brien (1989), 
and leads consequently to a lower limit of the calculated 
i^F yields. 

3. Predicted ^^F yields from individual WR stars 

As discussed by Arnould et al. (1999) on groimds of the 
NACRE rates, ^^F could be overproduced (with respect 
to solar) by the CNO cycle only at temperatures around 
15 X 10^ K, the exact level of this overproduction remain- 
ing poorly predictable, however, in view of remaining rate 
uncertainties. This conclusion contradicts the one derived 
from the use of the rates recommended by Caughlan & 
Fowler (1988), in which case fiuorine can never emerge in 
significant amounts from the CNO burning. As the latter 
rates are adopted in our calculations, the CNO zones of 
the computed model stars are depleted in ^'-^F. This trans- 
lates directly into a decrease of the ^^F mass fraction Xf g 
at the stellar surfaces when the ^^F-dcplcted CNO ashes 
are uncovered by mass loss (with the choice of the ordinate 
scales, the changes of fluorine abundance at the center and 
at the surface during the H-burning phase are not visible 
on Fig. 1). With the NACRE rates, it is expected that 
more would be present at the surface. However, it is 
also likely that this change is not able to affect drastically 
the predicted final yields, as these are dominated by the 
i^F made during the He-burning phase. 

In fact, as seen in Fig. 1, fiuorine builds up through 
i^N(a,7)i8F(/3+)i80(p,a)i5N(a,7)i9F during the early 
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Fig. 1. Evolution of the total mass Mtot) of the mass of the convective core Mconv, and of the central (^ig) and surface 
(Xfg) ^"F mass fractions for the 60 M0 model stars with mctallicities Z — 0.008, 0.020 and 0.040 during the end of 
the H-burning stage and the whole He-burning phase. The initial ^^F mass fraction is assumed to relate to the solar 
value XjI by X°g(Z) = (ZIZQ)Xf^. The spectroscopic types encountered during the evolution are indicated on the 
right of the figure: OV for 0-type main sequence stars, LBV for Luminous Blue Variables, WNL, WNE and WC for 
the different classes of WR stars. Note the different ordinate scales on the left and on the right of the figure. 
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phase of core He-burning. However, at the end of He- 
burning, ^^F (a , p) ^^Ne is responsible for a significant 
^^F destruction. Thus, material experiencing the whole 
He-burning episode cannot be -"^^F-enriched. In contrast, 
in massive stars going through the WR stage (initial mass 
Mi ^ 25 M0 for Z = 0.02, M; ^ 35 Mq for Z = 0.008; 



see Maeder & Meynet 1994), 
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F synthesized early 



during the core He-burning phase is ejected into the in- 
terstellar medium by stellar winds before its destruction. 
Indeed, Fig. 1 exhibits an increase of Xfg when the He- 
burning products appear at the surface during the WC 
phase. As a result, the ratio (Xfg(WC))/X^ of the aver- 
age -'^^F surface mass fraction during the whole WC phase 
to the solar system ^^F mass fraction takes values as high 
as about 55, 95 and 60 in the case of the 60 model 
stars with Z = 0.008, 0.02 and 0.04, respectively 

Figure 2 shows the ^^F "wind" yields for the computed 
stars (Mi, Z) with initial mass Mi and metallicity Z. These 
yields, noted pYg"^{Mi, Z), are equal to 



p™<i(Mi,Z) 



T(Mi,Z) 



M{MuZ,t)[Xl^{MuZ,t) 



■X%{Z)]dt, 



where T{Mi,Z) is the total lifetime of the star (Mi, Z), 
M{Mi, Z,t) its mass loss rate at age t, Xlg{Mi, Z,t) its 
-•^^F surface mass fraction at age t, and Xig{Z) its initial 
^^F mass fraction, assumed to relate to X{| by Xig{Z) = 
{Z/ZQ)X®g. These yields may be negative in case most of 
the ejected material has been depleted in fluorine. 

Figure 2 demonstrates that the highest yields are ob- 
tained for stars with Z = Zq and 40 <, Mi <, 85 M0. At 
lower mctallicities, the winds are indeed weaker, and thus 
uncover the He-burning core only for the most massive 
stars and when the ^^F has already been burnt. On the 
other hand, at higher mctallicities and for Mj > 85 M0, 
the H-burning core mass decreases so rapidly during the 
main sequence as a consequence of very strong stellar 
winds that the He-burning core becomes too small for be- 
ing uncovered by the stellar winds. 

The above discussion shows that the most important 
physical ingredient influencing the WR '^^F yields is the 
metallicity-dependent mass loss rates, quantities like con- 
vective core masses being less crucial in this respect. As 
a numerical example, the value for (Xfg(WC)}/Xfg rises 
from about 18 in the 60 Mq low mass loss rate model of 
Meynet & Arnould (1993a) to about 95 in the same model 
star computed in this paper with an increased M value. 
This high sensitivity to M might cast doubts on the relia- 
bility of the predicted ^^F yields. In fact, some confidence 
in the results presented in this paper may be gained by 
noting that our present choice of the mass loss rates al- 
lows to account for the variation with metallicity of the 
number ratio of WR to 0-type stars in regions of constant 
star formation rate (Maeder & Meynet 1994). 




40 60 80 100 
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(1) Fig. 2. Mass of ^''F ejected by the stellar winds (pf'J"'^ in 
Eq. 1) as a function of the initial mass and metallicity. 



4. Estimate of the contribution of WR stars to 
the galactic fluorine 

In order to evaluate the level of ^^F contamination by the 

winds of WR stars on a galactic scale, we use the p™'^ 
yields [Eq. (1)] in a very simple model of galactic chemical 
evolution making use of the closed box and instantaneous 
recycling approximations. We also suppose that only WR 
stars are able to affect the galactic ^^F budget through 
their winds, all other possible production or destruction 
sites being neglected. 

In such conditions, the ^^F mass fraction XiQ{t) in the 
galactic gas at time t is equal to (e.g. Tinsley 1980) 



Xi9{t) = yi9ln[l/a{t)], 



(2) 



where a{t) is the mass fraction of the gas in the Galaxy at 
time t, and yig is a representative time-independent ap- 
proximation of the net yield of a stellar generation defined 
by 



2/19 (i) 



I /-Ma 
i - R J Ml 



pYr{MuZ{t))^{Mi)dMu (3) 



where R is the "returned fraction" , Mi and M2 the lowest 
and highest mass of the stars going through the WR phase, 
and <i>(Afi) the initial mass function (IMF). It has to be 
noted that Eq. (2) would break down if the true time- 
(or Z-)dependent yi^{t) yields were used instead of yig. 
In order to evaluate the latter quantity, we notice that 
2/i9(Z(i)) values of about lO"'^, 5.7 10"'^ and 2.2 10"^ are 
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obtained for Z = 0.008, 0.02 and 0.04 if use is made of 
the "p™"^ values reported in Sect. 3 and of the (properly 
normalized) IMF derived by Kroupa et al. (1993). On such 
grounds, we just adopt the rough estimate j/ig « 3 10~^. 

If this approximation is used in conjunction with the 
value fj w 0.2 considered to characterize the solar neigh- 
bourhood at the time of the solar system formation 4.5 
billion years ago (see Prantzos & Aubert 1995, and refer- 
ences therein), Eq. (2) leads to w 510"'' in the local 
Z = Zq interstellar medium (to be compared with the 
solar system abundance of 4 10"''). Thus, our simple es- 
timate predicts that WR stars might account for most of 
the solar system ^^F content. Even larger ^^F quantities 
would be predicted with the use of the ^^N (n , p) ^"^C rate 
of Koehler & O'Brien (1989)! After having faced for long 
the problem of the underproduction of "'^^F, the theory 
of nucleosynthesis might now live with the danger of its 
predicted overabundance. If this is confirmed by further 
studies, constraints will obviously have to be put on one 
model or another. 

5. Implications of ^^F detection at high redshift 

Any ^^F present at high redshifts has to have been synthe- 
sized in massive stars only. Timmcs ct al. (1997) have ar- 
gued further that its detection at redshifts z > 1.5 would 
in fact be a signature of the z/-process in massive star 
explosions. The possibility of '^''F production by non- 
exploding WR stars might in fact weaken this statement, 
and blur the picture substantially. 

Of course, one has to acknowledge that the contribu- 
tion from WR stars at high redshifts may be reduced as 
a direct result of the lower metallicities that appear to 
characterize such regions. According to observations of 
Damped Lyman a systems (Pettini et al. 1997), the metal- 
licity at redshifts between 1.5 and 2 indeed lies around 
O.IZq. Such a reduced metallicity lowers the WR "'^^F 
yields for two reasons. First, the number of WR stars pre- 
dicted by non-rotating single star models is considerably 
reduced as a result of lower mass losses (Maedcr & Meynet 
1994). Second, the abundances of the CNO seeds that are 
needed for the secondary WR production are reduced 
as well. 

Even so, it would certainly be premature at this point 
to completely forget about the role of WR stars in a pos- 
sible enrichment of high-z material with "'^^F, and to relate 
it strictly with the j^-process.This is even more true as 
the predictions reported in this paper are based on single, 
non- rotating stellar models only. How binarity and/or ro- 
tation would change these results remains to be checked. 
At present, the published rotating evolutionary models 
leading to WR stars (Fliegner & Langer 1994, Meynet 
1998, 1999) make no predictions concerning the synthesis 
of fluorine. However they show that rotation favours an 
early entrance into the WR phase for a given mass, and de- 
creases the minimum initial mass for a star to go through 



a WR phase at a given metallicity. Moreover, the mixing 
induced by rotation opens up new nucleosynthetic chan- 
nels (see Heger 1998) whose importance for the scenario 
of fluorine production presented in this paper remains to 
be quantitatively assessed. Finally, let us note that the 
effects of rotation might be more important at low Z if, 
as suggested by Macder et al. (1999), the average rota- 
tion is faster at low metallicities. In such conditions, and 
in absence of quantitative calculations, one has to remain 
alert to the possibility of a signiflcant contamination of low 
metallicity high redshift regions by the ^^F-loaded wind of 
WR stars. 

Clearly, observations of -'^^F at high redshift, if possible 
at all, would be decisive in order to answer the question of 
the very production mechanism of this element. An impor- 
tant distinguishing feature would be the primary nature 
of the observed ^^F, as predicted by the j^-process, or its 
secondary behaviour, as expected from the thermonuclear 
model discussed in this paper. 

6. Conclusion 

Detailed stellar model predictions made in the framework 
of a very rough model for the chemical evolution of the 
solar neighbourhood leads to the conclusion that non- 
exploding non-rotating single WR stars alone could ac- 
count for the solar ^^F content. This conclusion remains 
to be ascertained by the adoption of a more realistic galac- 
tic evolution model. Still, it appears likely that the con- 
sidered WR stars might be signiflcant, and even possi- 
bly dominant, galactic ^^F contributors. In addition, they 
might well be responsible for a ^^F enrichment, if any, 
of high-rcdshift (z > 1.5) low-metallicity regions (~ 0.1 
Zq). Further predictions are eagerly awaited for rotating, 
as well as binary, WR stars. 

Finally, let us stress that the most direct test of the 
present model would be the measurement of the abun- 
dance of fluorine in the wind of WC stars. It remains to 
be seen if such observations are really feasible. 
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